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Abstract: In this paper we present a variety of nonlinear controllers for the magnetic bearing that ensure both 
stability and robustness. We utilize techniques of discontinuous control to design novel control laws for the magnetic 
bearing. We present in particular sliding mode controllers, time optimal controllers, winding algorithm based con- 
trollers, nested switching controllers, fractional controllers and synchronous switching controllers for the magnetic 
bearing. We show existence of solutions to systems governed by discontinuous control laws, and prove stability and 
robustness of the chosen control laws in a rigorous setting. We design sliding mode observers for the magnetic bearing, 
and prove the convergence of the state estimates to their true values. We present simulation results of the performance 
of the magnetic bearing subject to the aforementioned control laws, and conclude with comments on design. 


1 Introduction 


Magnetic bearings require high precision control systems to ensure adequate stability, stiffness and robustness. Mag- 
netic bearing controllers are typically linear controllers that perform pole placement through state feedback. Nonlinear 
control techniques on the other hand provide surprisingly simple, yet stable and robust controllers. It is our con- 
tention that nonlinearities can be introduced and exploited in a controlled manner to provide significant performance 
enhancements without increasing the complexity of control. To this end we present the following control techniques 
and controllers for the magnetic bearing. 

• Sliding mode control of magnetic bearings. 

• Minimum time control of magnetic bearings. 

• Winding algorithm based control of magnetic bearings. 

• Nested switching control of magnetic bearings. 

• Synchronous control of magnetic bearings. 

• Fractional control of magnetic bearings. 

• Sliding mode observers for magnetic bearings. 

Each of the aforementioned control methodologies utilizes discontinuous control of the magnetic bearing. We 
show rigorously the existence of solutions to these highly nonlinear systems and prove the stability and convergence 
of the magnetic bearing system subject to each of these control laws. 

The organization of this paper is as follows. In the first section we present the basics of the theory of differential 
equations with discontinuous righthand sides. The second section presents the dynamical equations of the magnetic 
bearing. The third section presents the theory, proof and simulation of a sliding mode controller of a magnetic 
bearing. In the fourth section we present the theory, proof and simulation of a minimum time controller of a 
magnetic bearing. Section V of this paper presents winding algorithm based control of the magnetic bearing. The 
sixth section presents the theory, proof and simulation of a nested switching controller of a magnetic bearing and is 
followed by the synchronous controller in the eighth section. The ninth section of this paper presents a fractional 
controller and a conjecture pertaining to the stability of the fractional controller. We conclude this paper with a 
presentation of sliding mode observers in the tenth section. The appendix lists facts and definitions from real analysis 
useful towards understanding the mathematics presented in this paper. 
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2 Differential Equations with Discontinuous Righthand Sides 

As a prelude we compare and classify ordinary differential equations based on the nature of their right hand sides 
Consider a differential equation of the following form. 


i = /(*,t) 

x(0) = io 

x € 3?" t 6 5R+ 
/(x, f) : S?" x 5R+ 


■ 5 ?" 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


The smoothness assumptions on f(x,t) determine the kind of differential system referred to by (1). 

The three major kinds of differential systems are 

1. Cauchy Differential Systems: In the domain D of the ( x,t ) space, 

• /(x,t) : 9?" x 9?+ — ♦ is continuous in x. 

• /(x, 0 : is continuous in i. 

2. Caratheodory Differential Systems: In the domain D of the (i,<) space, 

• /(x,*) : 5R n x 9£ + — ► ?R n is continuous in x. 

• /(x,t) : 5R n x is discontinuous in t on sets of zero measure. 

3. Filippov Differential Systems: In the domain D of the (x,t) space, 

• /(x,f) : 9R n x is discontinuous in x and i on sets of zero measure. 

The nature of the righthand sides indicates the kinds of solutions (strong or weak) that exist for the d.fferential 
system. From a control systems engineering standpoint, the use of discontinuous controls necess.tates the need fora 
careful stability analysis owing to the nature of the solutions that exist for such systems. 

2.1 Filippov Differential Systems 

In this section we will develop solution concepts and conditions for existence of solutions to 

with discontinuous right hand sides. Such equations represent physical systems governed by switching behaviours. 

Instead of describing solutions for differential equations with discontinuous right hand sides, we will consider 
differential inclusions which include the said discontinuity as a special case. We will then describe gen«a.hz 
solution concepts for these differential inclusions, and will present conditions for ex.stense of generalized solutions 

differential inclusions. 

We consider Filippov Differential Systems of the following form. 

(6) 


i = /(*,<) 

x(< = 0) = xo 

X e 3?" t € »+ 

f(x,t) : x »+ — 9f" 


( 7 ) 

( 8 ) 
(9) 

( 10 ) 


where in the domain D of the (x,t) space, 

• f(x,t) : 9J n X 5?+ — 1 • 5?" is discontinuous in x 6 D on sets of zero measure. 

• /(x l) : 9i" x Si-t- — * 5? n is discontinuous in t € D on sets of zero measure. 

• f(x,t) : 5R n x 5R+ — ► SR" is measurable in t € D for each x € D. 

• ||/(x, t)|| < Kj{t) V (x,t) € D where Kj(t) : SR+ — SR is summable. 

The aforementioned conditions on the function /(«.*)*" * - *" are also Ca,,ed FUipP ° V conditi ° nS ' 

We will now consider a differential inclusion that adequately describes the discontinuous system^ Though 
r . ftr ll • x — SR" of equation (6) is undefined on sets of zero measure, we choose instead to represent 
r / n ffx tl »" x SR + by a set valued map on such sets of zero measure. That is to say, if for 

the function f(x ) . + , * t *x £ « we f orma ll y define the function to be set valued 

instance the function is undefined at a point ^x , t ; fc * x •*+, c J 

,, „: nt ir* Indeed depending on the set-value attributed to the function at the point (x , ), 
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discontinuity. To justify the use of the inclusion, we must show that given any arbitrary e £ there exists a small 
enough 6 € neighbourhood of the point of discontinuity, such that, the trajectories of the differential equation 
in this 6 neighbourhood are e close to the solutions of the differential inclusion. Furthermore, as the size of the set 
containing the point of discontinuity shrinks to zero, that is 6 — 0, the solutions of the differential equation tend 
to the solution of the differential inclusion. That is to say, that the trajectories of the differential equation weakly 
converge to the solution of the differential inclusion. We will say more about this later. 

Indeed, given a discontinuous differential system of the form (6), henceforth we will replace it (whenever 
possible) with a differential inclusion of the following form. 


X £ F(X y t) 

x(t = 0) = Xo 

x 6 t e 

F{x,t) : 9T X 9?+ — S £ 

where 5 is a set in 9?" and in the domain D of the (x,f) space, 

• the set valued map F(x,t) : 9? n x 9?+ — ► 5 £ is tipper semi- continuous. 

• The Range[T{x,t)] £ 5?" is compact and convex. 


(11) 

( 12 ) 

(13) 

(14) 

(15) 


Comment 2.1 The definition of the inclusion T{x,t) is such that it is single-valued in the domain of continuity of 
the function f(x,t); indeed it is equal to f(x,t) in the domains of continuity , but is set valued in the domains of 
discontinuity of /(x,<). 


Comment 2.2 It is important to note the properties of the set S £ 9*? n which will be used for the existence of solutions. 
We now formally define the solution of a Filippov differential system. 

Filippov Solution Concept: An absolutely continuous vector function s(t) : 9£+ — is defined to be a 
Filippov solution of the Ftlippov differential system (11) if for almost all t £ D } 

where 


in the domains of continuity ( 17 ) 

*■(«(**). o = n n convex-hul\B{x, 6) — N, t) (lg) 

< 5 > 0 ,jN = 0 

and f\jv=o d enotes the intersection over all sets N of Lebesgue measure zero where the function f{x y t) is either 
undefined or discontinuous. 

Comment 2.3 In the domains of continuity of f{x,t) : 9?+ — 9? n , the inclusion F(x y t) is the same as the function 
and therefore the set operation £ in equation (158) must be replaced with the strict equality = 

The utility of the Filippov solution concept is that it is indeed the limit of solutions to (6) averaged over neigh- 
bourhoods of diminishing size. The key point to be understood is that the Filippov trajectories of the discontinuous 
system remain close to the true trajectories. 

As is evidenced in the proofs of the Cauchy and Caratheodory systems, the method of constructing solutions 
to differential equations begins by constructing sequences of approximating solutions, and then ensuring that the 
approximations converge in some sense. 

We now state the theorem that guarantees the local existence of Filippov solutions. 

Theorem 2.1 Local Existence Of Filippov Solutions To Filippov Differential Systems 
Given (Gl) A Fi/ippot> differential system of the form (11). 

If (II) The domain D where }(x y i) is specified for almost all t 
(x,t) £ 9£ n x 3£+ : | \x — Xo|| < K x and t < K t 

(12) f(x,t) : 5? n x 9?+ — ► is measurable in t £ 9?+ for all x £ 9?" 

(13) ||/(x,/)|| < hf{t) V (x t t) £ D where K /(/) : 9?+ — ► 9? is summable. Furthermore there exists Kj £ 9? + 
such that Kj > |A'/(<)| V i £ D 

(14) The differential inclusion in (11), T(x,t) : 9?" x 5R+ — 5 € 3T, where S is a set in 9? n and in the domain 
D of the (x,t) space satisfies the following two assumptions. 

* the set valued map F(x y t) : 9£ n x — ► S £ 9? n is upper semi-continuous. 

* the set S £ is compact and convex. 

Then (Tl) The differential system (11) has at least one Filippov solution s(t) : — ► 9? n for t < min(K t , —*■) 

satisfying the initial condition s(0) = xq. 
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3 Dynamic Equations of the Magnetic Bearing 


The dynamic equations of the magnetic bearing [RG93] may be written as follows. 


xi _ 

X2 


04X4 
O 4 x4 


h x4 


= [^4x4 04 x4 ] 


Xl 

X2 

*1 

A 


+ 


0 

B 


where, 


x\ € & 
x' 2 € ^ 


u € 


y e® A 


A : ft* — ft 4xi 


i e ft 4 * 4 


B € ft 4 * 4 


*1 

x? 

*1 

*2 
*2 
*2 
L x 2 


L ST J 


’ 0 0 
0 0 
0 0 
. 0 0 

" 1 0 
0 1 
0 0 
0 0 


0 0 

0 0 

0 — aw 
aw 0 

0 0 ' 

0 0 

1 0 
0 1 


ri i 0 0 

mm 

0 0 ^-- 
mm 

0 0 — /l h 

f, — f 2 0 0 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


The state variables xi € 5ft 4 are the generalized positions and rotational angles while the state variables 12 € 5ft are 
the generalized linear and angular velocities. 

Now choose a control law of the following form. 


u = 


v = 


B 1 [-A[u)]x 2 + v] 


(29) 

(30) 


where the control inputs t. € 5ft 4 will be specified later. Substituting control law (29) in the dynamical equations (19) 
- (20), we arrive at the decoupled form of the state equations written as follows. 


i'l = *7 

i' 2 = v' 1 = 1,2 4 


(31) 

(32) 


4 Sliding Mode Control of the Magnetic Bearing 

In this section, we specialize the theory of discontinuous systems to a special class of systems of the following form. 
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X 


= /+(x) for [x : s(x) > 0] (33) 

= /-(s) for [x : s(x) < 0] (34) 

where x £ and f(x) : — ► 5?” and s(x) : 5?” — ► -ft. Note that S = {x : s(x) — 0} is a manifold of dimension 

n — 1. This manifold 5 is called the sliding manifold or sliding surface. The dynamics of the system on this manifold 
S is called the sliding dynamics or sliding modes of the system. The design of the manifold S is such that it is globally 
attractive, and trajectories commencing from arbitrary initial conditions reach S in finite time. Furthermore, the 
dynamics on 5 achieves the control objective. 

Local existence of solutions is verified by modelling the system represented by equations (33) - (34) by the appro- 
priate differential inclusions and verifying whether the inclusion satisfies the hypotheses of the theorem concerning 
local existence of Filippov solutions. 

Uniqueness, in the sense of the Filippov solution is shown if either (x) < 0 or f-( x ) > 0- This is 

shown in [SS83], [Fil88], [Fil6 1 ] . The physical interpretation of these conditions is simply that the trajectories of the 


system are always directed towards S y thus rendering it attractive. 

Example 4.1 

x = —ksgn[x] (35) 

sgn[x] = 1 ifx > 0 (36) 

sgn[x) — — 1 ifx < 0 (37) 

Modelling the system (35) by a simple differential inclusion, we rewrite (35) as 

i £ T(x) (38) 

where 

T{x) = sgn[x] ifx ^ 0 (39) 

T{x) £ [—1,1] i/ x = 0 (40) 


The inclusion in (38) is closed, bounded, convex and upper semicontinuous and therefore by the theorem on existence 
of Filippov solutions, Filippov solutions exist for this system. 

The sliding modes of a system, defined to be the Filippov solutions to the system on the manifold 5, are calculated 
by performing Filippov averaging, which is a convex combination of dynamics on either side of the manifold 5. Indeed, 
by dynamics on either side of the manifold 5, we merely refer to /+(x) and /_(x). The simple extension of the notion 
of sliding manifolds to non-autonomous systems is shown in [SS83]. 

While the theory of existense of solutions has been developed for general nonlinear systems with discontinuous 
controls, the methodology to design sliding mode controls to achieve stabilization or tracking is well understood 
only for a restricted class of systems [SRS91]. In the following sections, we will present the theory for Linear Time 
Invariant Systems - SISO and MIMO. 


4.1 Sliding Mode Design For LTI Systems 

Consider linear time invariant systems represented by the following equations 

x = Ax + Bu (41) 

where x £ A 6 9? nxn , B 6 9? mxn and the controls u £ We will now prescribe the sliding mode controller 

design procedure in a sequence of steps. 

Step 1. 

Check to see if the system is completely controllable. If the system is not completely controllable, a sliding mode 
controller cannot be designed. 

Step 2. 

If the system is completely controllable, find a linear transformation of the state that recasts the system in the 
controllable canonical form. That is find a transformation 

x = Tx T £ 5i nXn (42) 

such that the state equations are of the form 
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Step 3. 

We define 5(x) : — 3? as 


5(x) = ajXi + « 2 X 2 + * — h a n _ix n _i + x n ' ' 

where the coefficients a, « = 1, 2, . . . , n - 1 of (44) are such that the polynomial S(x) is a Hurwitz polynomial. 
Furthermore, note that S = 0 is an n - 1 dimensional manifold, called the sliding surface. 

Indeed now choose the control input u to be 


ll(t) = — 61 Xi — 62X2 • ‘ — fenXn “ 

v(t) = — ai X2 — a>2 X3 — * * * — a n— 1 — fcs<jr7i[5(x)] 


(45) 

(46) 


Choice of control u enables us to rewrite the system in the form 


Xi = x 2 
X2 = 

X n _j = — Ui X 1 — 02^2 — • • * — flyj — 1 Xn — 1 “H ‘^(• c ) 

S(x) = -fc$ 0 n[s(x)] 

It is easy to show that Filippov solutions exist, and that S(x) = 0 is reached in finite time from arbitrary initial 
conditions. Furthermore on the n — 1 dimensional manifold 5 = 0, the reduced order dynamics is exponentially 
stable. Consequently global exponential stability of the system is shown. 

The choice of discontinuous input induces chatter in the system. To reduce the chatter, we utilize various 
regularizations and smoothings of the discontinuous sgn function. The common smoothing technique is the use of 
the saturation function, which is presented in [SS83]. 

We now present a choice of continuous control input that enables us to reach the sliding surface S = 0 in finite 
time. Indeed, consider the control given by 


(47) 

(48) 

(49) 

(50) 


u(t) = — X l — 6212 ’ ‘ Cj ( t ) 

v(t) = - 01 X 2 - 02*3 a„-i*„ - *|S(x)|"Sffn[S(x)] 

m > 1 


(51) 

(52) 

(53) 


Such a choice of control u enables us to recast the system equations in the form 


X] = *2 


*2 


Xn — 1 

S(x) 


*3 

— Ol X J — 02 X 2 0 ,,-lIn-l + S(x) 

— *l 5 (*)| ” *!?»[*(*)] 


(54) 

(55) 

(56) 

(57) 


It is easy to show that Filippov solutions exist, and that the n — 1 dimensional manifold S(x) — 0 is reached in 
finite time. Furthermore on the n - 1 dimensional manifold given by 5 = 0, we see that the reduced order dynamics 
is exponentially stable. Consequently global exponential stability of the system is shown. Phis control law u is 
interesting in that it is continuous, but not differentiable. 


Comment 4.1 The disturbance rejection properties of the discontinuous control law are significantly better than that 
of the continuous control law . This indeed is the design tradeoff involved in designing continuous control laws . 

Comment 4.2 The extension of the sliding mode control techniques to controllable MIMO systems that are decouplable 
is trivial. Once the system equations are transformed into decoupled systems , each of which is in the controllable 
canonical form, we apply the design method outlined earlier to design sliding surfaces for the decoupled system. Note 
however that sliding occurs not at the individual surfaces , but at the intersection of all these surfaces. 


Theorem 4.1 Sliding mode control of a magnetic bearing system. 

Given (Gl) A magnetic bearing system of the form (31) - (32). 

If (II) The controls v' i = 1,2, ... ,4 are chosen as 

v‘*l,ding = -“ 1*2 - k * + * 2 ] ( 58 

Then (Tl ) Filippov solutions exist for the system (31) ■ (32) subject to the control law (58). 

(T2) The trajectories of the system (31) - (32) subject to the control law (58) reach the ortgtn in fimte time. 
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Figure 1: Sliding Mode Control of a Single Axis of a Magnetic Bearing 

Proof: 4 t> To show the existence of generalized Filippov solutions we first note that the dynamical system 
represented by equations (31) - (32) subject to the control law (58) can be modelled by the following inclusion. 


i\ = -aJzJ+S (59) 

S e T\x) (60) 

where the inclusions «F'(x) : 3ft — ► [— k,k] are specified as 

T\x) = — k * $gn[a\x\ + x* 2 ] if ||x ||2 > 0 (61) 

£ [— fc, fc] else (62) 

i = 1,2,..., 4 (63) 

The inclusions T x (x) i = 1 , 2, . . . , 4 are 


• closed, bounded, convex and uppersemicontinuous. 

Invoking the theorem on the existence of generalized Filippov solutions, we conclude that Filippov solutions exist for 
the system (31) - (32) subject to the control law (58). 

Stability and robustness of the magnetic bearing follow from the earlier discussions. 

The phase portrait of trajectories subject to the sliding mode control is given below. 


5 Minimum Time Control of the Magnetic Bearing 


It is many times desirable in a magnetic bearing to choose a control law to perform regulation in minimum time. Such 
minimum time regulation ensures good response to impulsive perturbation forces. To achieve regulation in minimum 
time, we formulate the optimal control problem as specified in [AEB75]. 

Consider the minimum time optimal control problem with the functional to be minimized, given by 

J = / dt (64) 

Jo 


Theorem 
Given (Gt) 

(G*) 
If (II) 


5.1 Minimum time control of a magnetic bearing system. 
A magnetic bearing system of the form (31) - (32). 

A functional to be extremized of the form (64). 

The controls v* i — 1, 2, . . . , 4 are chosen as 


Vopttmal 


sgn[x\ + if\x\ + ^\>0 

-sgn[x' 2 ] if\x\ + = 0 


(65) 
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Then (Tl) The trajectories of the system (Si) - (32) subject to the control law (65) reach the origin tn minimum 
time. 


Proof: * t> Using standard methods of optimal control [AEB75], we write down the Hamiltonian function 
H(x t u, A, <) as 

H(X, U, A, f) = 1 + Al x\ + A 2 v' 

Inspection of equation (66) reveals that the control that minimizes the Hamiltonian is given by 
v x = — s0n[A2]vm OI 

where t4 OI is the maximum permissible value of control. Without loss of generality, we will assume that v' max = 1. 
where Ai and A 2 are the co-state variables. The co-state equations are given by 


Ai = 0 

A 2 = —Ai 


( 68 ) 

(69) 


Integrating the co-state equations yields 
A 2 (t) = -Ai(0)t- A 2 (0) 

Therefore the optimal control is given as 
v l = sgn[— Ai(0)< — A 2 (0)] 

The control can assume only two values +1 or — 1. When v* = +1, we integrate the state equations to obtain 


x\ (t) 
x\(t) 


t + £ 2 ( 0 ) 

— -f X2(0)< 4- z*i(0) 
Eliminating t we obtain 



+ * 1 ( 0 ) 


[4] 2 W 

2 


Similarly, when v' = -1, integrating the state equations we obtain 
x\ (t) = -t + *i(0) 

x\(t) = — — + + *1(0) 

Eliminating / we obtain 


x y = - 


Ltt +rfw .acia 


These curves describe a family of parabolas, whose switching curve may be written as 

x' 2 \x' 2 \ 


S(x,t) = x\ + 


(70) 

(71) 

(72) 

(73) 

(74) 

(75) 

(76) 

(77) 

(78) 

(79) 


(80) 


In terms of the switching curve, the control v’ optima i may be written as 


_ f -sgn[x[ + 
V ° pUrna ' \ -sgnte) 


if \x\ + 

if 1*1 + 


2 1 

ilfail — 
2 


> 0 

= 0 


(81) 


^The phase portrait of trajectories subject to the optimal control v‘ optttnal is given below. Note the trajectories 
converging to the switching curve, which is nonlinear (while the switching curve in conventional sliding mode systems 
is linear). The chosen control gains are 


it, = 1 

k 2 = 2 


(82) 

(83) 
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Figure 2: Minimum Time Control of a Single Axis of a Magnetic Bearing 


6 Winding Algorithm Based Control of the Magnetic Bearing 

The winding algorithm was introduced by [Kor68], [Pra92] and makes use of continuous switching between the surfaces 
ij = 0 and x 2 — 0 to reach the origin. The interesting feature of this control technique is that the control has two 
switches. One switch is used to change the direction, and the other is used to change the magnitude. By repeatedly 
switching between the surfaces x\ = 0 and x 2 = 0, we wind closer to the origin. 


Theorem 6.1 Magnetic bearing control utilizing the winding algorithm. 
Given (Gl) A magnetic bearing system of the form (31) - (32). 

If (ll) The controls v‘ i = 1,2, ... ,4 are chosen as 


Then (Tl) 

(T2) 


V winding — k 2 sgn[x 2 ] k\ > k 2 > 0 (^4) 

Filippov solutions exist for the system (31) - (32) subject to the control law (84). 

The trajectories of the system (31) - (32) subject to the control law (84) wind to the origin in finite time. 


Proof: 4 t> To show the existence of generalized Filippov solutions we first note that the dynamical system 
represented by equations (31) - (32) subject to the control law (84) can be modelled by the following inclusion. 

x\ = x 2 (85) 

6 ( 86 ) 

where the inclusions T'(x) : 5R -» [-(* a + k 2 ), (*] + fc 2 )] are specified as 


T'(x) = -fcjsj/nfxJ] - fc 2 sgn[z 2 ] if ||z|| 2 > 0 

€ [— (fci + k 2 ), (fci + £ 2 )] else 

* = 1, 2, ... ,4 

The inclusions -F'(x) j = 1, 2, ... ,4 are 


(87) 

( 88 ) 
(89) 


• closed, bounded, convex and uppersemicontinuous. 

Invoking the theorem on the existence of generalized Filippov solutions, we conclude that Filippov solutions exist for 
the system (31) - (32) subject to the control law (84). 

Let us first prove the stability and finite time stabilization of the algorithm. To show stability, we use the extended 
Lyapunov theorem, [AKP91] proofs for which may be found in [AC84]. The theorem is primarily used to conclude 
weak- stability oi differential inclusions by investigating generalized gradients of non-differentiable Lyapunov functions. 
A brief statement of the theorem would be as follows. 

Given a differential inclusion x e F(x,t) and a nondifferentiable Lyapunov function V(x). If for every element 
v in the generalized gradient of V , there exists at least one element f £ F(x,t), such that L F V < 0, then the zero- 
solution is weakly asymptotically stable.. Indeed, weak asymptotic stability is the best we could hope for when dealing 
with set-valued differential inclusions. 
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Now consider the system (31) - (32) subject to the controls u w ,„ding- The system equations are 
i 1 = x 2 

x 2 = — JfciSjfn[x‘i] — 

Consider a candidate Lyapunov function 


(90) 

(91) 


v= i* i|+ §r 

The derivative for xj.zj ^ 0 is given by 


V = - 

< 0 


k 2 x\ 


(92) 

(93) 

(94) 


Therefore xj - 0, and the reduced dynamics is such that xj - 0. However, when x) = 0, it is clear that we have 
to investigate the properties of the generalized gradient of V. However, it is obvious that when x , - 0, for every 
element t> of the generalized gradient of V, (which in this case happens to be any real number in (-1,1) ) there ews s 
an element of the inclusion F(x,t) (indeed, choose / = v) such that the generalized gradient of V along the flow 
of the inclusion F(x,t) is negative definite. The conditions of the generalized Lyapunov theorem are satisfied, and 

Finite time is shown by considering the state equations of the planar dynamical system in the various quadrants^ 
Indeed, if the portrait of the system were to be drawn with x) along the x axis and xj along the y axis, we wou 
note the following. 


x 


i 


i 

x 2 


■ [*z ] 2 

Jtl + k 2 

*Jfci -*z 


in the first and third quadrants 
in the second and fourth quadrants 


(95) 

(96) 


Every instance the trajectory moves from the first quadrant through the fourth quadrant to hit the y axis, we see a 
contraction occurring in the magnitude of £2 following manner. 


[li]2(<i)= ir+i7 [l * l2(0) t97) 

From the third quadrant through the second to strike the y axis again, we see the following contraction. 

(Xj] 2 (/2)=^^ [ ^ ]2 (fl) (98) 

The state trajectory therefore winds to the origin. <!♦ 

The phase portrait of the planar dynamical system subject to the winding algorithm is illustrated below. Note 
the very interesting way in which the state trajectories wind to the origin. The values of chosen control gains are 


*1 = 2 

k 2 = 1 


(99) 

( 100 ) 


7 Nested Switching Control of the Magnetic Bearing 

Nested switching controls work weU for planar dynamical systems [Pra92]. The basic approach is to permit chatter 
about the dual sliding surfaces x) = 0 and xj = 0. It is to be noted that chatter for multiple sliding surfaces is the 
equivalent of limit-cycle like behaviour. Consequently, by utilizing multiple sliding surfaces, and nondifferentiable 
controls, we are willing to tolerate limit-cycle like behaviour at the origin. Indeed, the problems associated with e hm- 
inating chatter in one-dimensional systems naturally extend to the higher order systems also. The use of saturation 
functions to perform nested switching is an extension of the idea of using saturation functions in one-dimensional 

systems, to many dimensions. 

Theorem 7.1 Magnetic bearing control utilizing the winding algorithm. 

Given (Gl) A magnetic bearing system of the form (31) - (32). 

If (II) The controls v' t=l,2,...,4 are chosen as 

t Wed = -k 2 sgn[x 2 - feisgufi)]] ( 101 ) 
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Winding Algorithm 



Figure 3; Winding Algorithm Based Control of a Single Axis of a Magnetic Bearing 


Then (Tl) Filippov solutions exist for the system (31) - (32) subject to the control law (101 ). 

(T2) The trajectories of the system (31) - (32) subject to the control law (101) reach the origin in finite time . 


Proof: 4 t> 

To show the existence of generalized Filippov solutions we first note that the dynamical system represented by 
equations (31) - (32) subject to the control law (101) can be modelled by the following inclusion. 


- *2 ( 102 ) 
x x 2 € F\x) (103) 

where the inclusions T'{x) : B? — [-(k x + * 2 ), (*i + *2)] are specified as 

T'{x) = -k 2 sgn[x l 2 - kisgn[x\]] if ||x|| 2 > 0 (104) 

C [— £2,^2] else (105) 

* = 1, 2, ... ,4 (106) 

The inclusions ^‘(x) i = 1, 2, . . . , 4 are 

• closed, bounded, convex and uppersemicontinuous. 


Invoking the theorem on the existence of generalized Filippov solutions, we conclude that Filippov solutions exist for 
the system (31) - (32) subject to the control law (84). 

Consider the system (31)-(32) subject to the nested switching control law given by 


x\ — 

*2 

(107) 

X 2 = 

— k2sgn[x % 2 T fcis^nfxj] 

(108) 

Now consider 

the following nondifferentiable Lyapunov function 


1/ 

[x* 2 -f Arisen]*!]] 2 


V — 

2 

(109) 

V = 

[x 2 + ki3gn[x\]][x 2 + 0]if \x\ \ > 0 

(110) 

= 

-k 2 \x x 2 -f fciS0n[x l 1 ]| 

(111) 

< 

0 

(112) 


Therefore x\ -*■ - kisgn[x\ ]. Indeed, it is easy to see that this happens in finite time. As in finite time x' 2 = 
-fciscjnfxJ]; now consider the Lyapunov function 


Vi 



(113) 
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Nested Switching 
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Figure 4: Nested Switching Control of a Single Axis of a Magnetic Bearing 


Vi = 

= xj[— JbiS9n[i‘i]] in finite time 

< 

< 0 


(114) 

(115) 

(116) 
(117) 


It is dear that *j - 0 in finite time. However, when x* a = 0, and is not equal to 0. I his is where 

chatter commences, and the system limit cycles between the surfaces x\ = 0 and x 2 - S “ ch h ”“ t c y c ‘"8 

behaviour is present as the gain ki is not slowly reduced as x*. - 0. Indeed if the mu tiphcand of sgn[z was to 
decrease in magnitude and finaUy equal 0 when x\ = 0, we can expect xi to also be equal to 0 without chatter. This 
indeed is the principle behind using saturation functions as opposed to sgn functions in nested control. We will now 
show an extension of this method, without using saturation functions. 

We now try to eliminate the problem of limit cycling between switching surfaces that was mentioned earlier. e 

do this using the switching control law mentioned earlier which is of the form. 


u, w ,tch,ug= -hsgn[x' 2 + fci |*1|™ S9«[*i]] 


(118) 


Denote S = x 2 + fci|*il m s S n t*i]- 
derivative of S may be written as 

S = -fc2S£fn[S] + k\ - 


Note that S is not differentiable at x\ = 0. However, almost everywhere, the 


li-dr 


(119) 


By choosing a large value of fe, we hope to swamp the term fci— Indeed, only in cases when this is possible, 
it is possible to conclude that 


x\ = -*i|x‘,| ™sjn[xi] 


( 120 ) 


And the conclusions of the previous section follow, without the limit cycle behaviour. 

The phase portrait shown below illustrates the properties of the control law. The values of chosen gains are 


k x = 0.5 

*2 = 5 

m = 2 


( 121 ) 

( 122 ) 

(123) 


For the same values of control gains, it is possible to choose a higher order fractional index, and the resulting 
phase portrait is shown below. 
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Figure 5: Nested Higher Order Switching Control of a Single Axis of Magnetic Bearing 

8 Synchronous Sliding Control of the Magnetic Bearing 

In this section, we present an interesting property of a modified vector sliding mode control law [AKP92a], [AKP92b], 
[AKP 5], and its possible application. The property of this modified vector sliding mode control law is such that 
it achieves simultaneous regulation for a group of n scalar systems with n inputs [SMSar], [SMS93], [SMS 7].This 
control technique has interesting implications for the magnetic bearing. Using this control technique it is possible to 
regulate the states of the multivarious magnetic bearings to the origin synchronously thus eliminating the possibility 
of inducing overshoot torques. The control law has the interesting property that it is a closed loop control law which 
can be prescribed without explicit reference to the initial conditions of the system. The law is interesting in that it 
introduces coupling between decoupled systems to achieve the synchronization objective. We present the basic theory 
of synchronous sliding control, and later specialize it to the case of the magnetic bearing. 


8,1 Synchronous Sliding 

Consider a group of n scalar decoupled systems of the form 


Xl 


til 

Xn 


ti n 

Xi(0) " 


X\Q 

Zn(0) 


X n0 


(124) 


(125) 


where the states x, G 9ft, the controls m G 3ft t = 1, 2, . . . , n the initial conditions x i0 G 3ft i = 1, 2, . . . , n. With minor 
abuse of notation, we create a new state vector x G 9ft n , where x — [ x\ • • ■ x n 

The control objective is to regulate the states from non-zero initial conditions to the origin, in finite time. That 
is, that there exist instants of time t * < oo G 9ft+ * = 1, 2, . . . , n such that the following is true. 


x,(0 = 0Vt>**i = l,2,...,n (126) 

We choose n sliding mode control laws of the following form to ensure achievement of the control objective. 

ii\xi\^0 i = l,2, ...,n (127) 


where k t G 9ft+ 


Comment 8.1 We note here that the controls u, t = 1,2,... are decoupled, in that u t is a function only of x t . 
Also note that the time taken by each state x t i = 1,2, ...,n to reach the origin is a function of its initial value 
x, (0) i = 1, 2, . . . , n and the control gains k t i = 1, 2, . . . , n. 
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Comment 8.2 Also note in equation (127) we did not specify the control law at |x,| = 0 i = 1,2, Indeed, 

u, = — lt,syn[i,] |x,| ^ 0. We do not specify the control at |x,| = 0. As the control is not specified only on sets 
of zero measure, it does not affect the existence of Filippov solutions shown by modelling the system by a differential 
inclusion . 


We now present some interesting properties of a modified sliding mode control law that deliberately introduces 
coupling between the decoupled systems. We present proof of existence of solutions, proof of stability, and proof of 
synchronous finite time convergence for the modified sliding mode control law. In order to do so, we formalize the 
notion of synchronous finite time convergence. 

Definition 8.1 A set ofnd Z+ variables x,(t) : R+ - R * = 1, 2, . . . , n are said to reach the origin synchronously 

commencing from nonzero initial conditions r,(0) / 0 i = 1, 2 n if there exists an instant of time t* < oo € 3?+ 

such that the following is true . 

Xi(t) ? OVfCt* 

Xi(t) = 0 Vt>f* 

i = 1,2,.. .,» 

That is to say, that the states with nonzero initial conditions (an assumption we make without loss of generality) 
are regulated to 0 at the same instant of time t*. There are many practical applications where such synchronous 
regulation is important. A typical application is a multifingered robot hand that grips an object. It is important to 
ensure that the fingers touch the object synchronously and thus cause force closure without imparting motion to the 
object. We will say more about this later. 

It is possible to ensure synchronous motion using a simple sliding mode feedback where the control gains are chosen 
with explicit dependence on initial conditions. Indeed, given the initial conditions exactly, we choose a decoupled 
control law that uses the values of initial conditions to derive control gains that guarantee synchronous reaching of 
the origin. For the sake of completeness we state the control law as follows. 


(128) 

(129) 

(130) 


Theorem 8.1 Synchronous regulation with explicit dependence on initial conditions. 
Given ( Gl ) A nonlinear system of the form (12)) - (125). 

(G2) A control law of the form (127) 

If (II) k, t = 1 , 2, . . . , ti are chosen such that 


MgH _ I x j(°)I 

k, kj 


1, 2, . . . , n j — 1,2 n 


Then (Tl) Filippov solutions exist for the system (124) * ( 1 25) subject to the control law (127). 
(T2) The surfaces x, = 0 t = 1, 2 , ... n are reached synchronously at a time t* = . 


(131) 


Proof: 4 t> The proof is quite straightforward and utilizes standard facts from sliding mode control theory. 

The existence of Filippov solutions is shown using the fact that the modelling differential inclusions F,(j) : 3? — ♦ 
[-1, 1] are closed, bounded, convex and uppersemicontinuous. Note that T,(x) : — [-1,1] are defined as follows 


t ,( x ) = -fc'-ri ^ l*«l ^ 0 

|x,| 

€ [-1, 1] if |x,| = 0 


(132) 

(133) 


2 

Stability is shown using the candidate Lyapunov function V{x) : 9i n -► »+ given by F(z) = Y.,=i whose 

derivative along the flow of (124) - (125) is given by V = -£" =] |x.|. Indeed V is negative definite proving global 

exponential stability of the origin. 

Finally, the time taken to reach the origin is given by t* = i = 1, 2, . . . , n. Now using the assumption that 

j = 1, 2, . . . ,« j = 1, 2, . . . , n } we see that t* = • ■ • = = **• 

This completes the proof of the theorem. <d^ 


Comment 8.3 The control law is inelegant to implement as it explicitly depends on the initial conditions. It would be 
desirable to develop a state feedback control law that would achieve the same objective, but one whose control gains 
do not explicitly depend on initial conditions. 


We now propose a state feedback control law that would ensure synchronous regulation. 

Theorem 8.2 Synchronous regulation with state feedback. 

Given (Gl) A nonlinear system of the form (124) - (125). 
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( 134 ) 


If (II) The controls u t i = 1,2, ... ,n in equations (124) - (125) are chosen to be 


= -k* 


if \\x\\ 2 > 0 * = 1,2, ...,n 


11*11* = (135) 

1=1 

(136) 

where k* 6 5R+ 

Then (Tl) Filippov solutions exist for the system ( 124 ) - (125) subject to the control law (134). 

(T 2 ) The surfaces x» = 0 t = 1, 2, . . . n are reached synchronously at a time t* = l i J W l h where | |x(0)| |a is the 
2 -norm of the vector of initial conditions, given by ||x(0)|| 2 = [£" = , x?(0)]i 

Proof: 4 t> We prove the theorem in three steps. First we show existence of generalized Filippov solutions to 
the system (124) - (125) subject to the control law (134). We then show attractivity of the origin when subject to 
the control law using a simple Lyapunov argument. Finally we show the achievement of synchronous regulation, by 
explicitly computing the times taken to reach the origin. We first make the following comments. 

Comment 8.4 It is interesting to compare the control laws given by equations (127) and (134). While the control 
specified by (127) decouples the system entirely , the control specified by (134) introduces a coupling between the 
through the 2 -norm of the state vector ||x|| 2 . Furthermore , note that the control gains k* remain the same for all 
Ui i = 1 , 2, . . . , n. 

Comment 8.5 The discontinuous control law (134) is not defined at the origin, the same way the function S£n[(-)] : 
^ ' ” [ — Iil] 15 not defined when (•) = 0. But also note that the control law specified by (134) * 3 bounded by k * . Indeed, 

as Tlffe - 1 * = i.2, . - - , Ui < **1 = i, 2 ,...,™. 

Step 1: Existence Of Filippov Solutions 

To show the existence of generalized Filippov solutions we model the system (124) - (125) subject to the control 
law (134) by the following differential inclusion. 


where the inclusions F t (x) : $ — ► [— k *, k*] are specified as 

T,(x) = if ||i|| 2 > 0 (138) 

€ [-**,**] if ||x|| 2 = 0 (139) 

* = 1 . 2 , ...,n ( 140 ) 

The inclusions F,(x) i = 1 , 2 , . . . , n are 

• closed, bounded, convex and uppersemicontinuous. 

Invoking the theorem on the existence of generalized Filippov solutions, We conclude that Filippov solutions exist for 
the system ( 124 ) - ( 125 ) subject to the control law ( 134 ). 

Step 2: Attractivity Of The Origin 

Consider a candidate Lyapunov function V(x) : 9? n — ► 5R+ given by 


Differentiating V along the flow of (124) - (125) subject to the control law (134), we find 

r 1 


= ~k* 


= -**||z||2 if ||xf| 2 / 0 
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Negative definiteness of V confirms the global exponential stability of the origin. 

Step 3: Synchronous Reaching 

From system (124) - (125) subject to the control law (134) the following is true for any i, j 



dx t Xi 

dxj x } 

V *, j < n i ^ j ||x|| 2 ^ 0 

Solving (148), we obtain explicit expressions for constraints on state trajectories as 

ii(<) = i,j<n i ^ j ||x||2 ^ 0 

x,(UJ 

Using (150) in (146), we recast (146) in the form 


= - k * 

= -k* 

= -V 
= 




X, 


ED-i*!]* 

x, 




fl 2 +T' n 

L 1 ! ^ Z_-/fc= 1, k*. 7*(o) 


= -k* 


= -it* 


X, = 




*<[! + 

S'(°) 

x,(0) 


r (0)|| 3 


* = 1,2, ... ,n 


The righthandside of (157) is a real constant, and therefore the solution of (157) is given by 


(146) 

(147) 

(148) 

(149) 

(150) 

(151) 

(152) 

(153) 

(154) 

(155) 

(156) 

(157) 


*' < ‘ , = _ ‘*iR°5K ,+ ’' (0) i = 1,2 * 

From (158), we obtain the time t* taken by x,(f) *=1,2,... 
initial conditions by setting the righthand side of (158) to 0. 


(158) 

n to reach the origin, starting from arbitrary nonzero 


0 

t* 


,* *-(°) 
11 ^( 0)112 
Mm ■ 


t * + x,(0) 

= 1,2, ...,n 


Synchronous convergence of state trajectories commencing from nonzero initial conditions is thus shown 
concludes the proof of the theorem. <!♦ 


(159) 

(160) 
This 


8.2 Design Of Tracking Control Laws 

The control laws that we have developed are discontinuous. As a prelude to presenting tracking control laws that 
involve discontinuities, let us analyze a simple linear pole- placement control law from another perspective. Consider 
a system represented as a chain of integrators of the form, 


Xl = X2 
in— 1 Xn 

Xn = U 


(161) 

(162) 

(163) 

(164) 
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where the state vector x 6 9^ and the control input u £ 9?. Given a desired smooth trajectory £id(<) : — ► 3? to 

be tracked by the state x x we present a tracking control law that uses successive derivative of desired trajectories. We 
define recursively, a set of desired trajectories for the states as 


Xxd{t) 


dx x - lyd (t) 

dt 


*i-l[Xi(0 - Xi,d(*)] 


2,3 


(165) 


While we are given a desired trajectory to be tracked by the state Xi(t), we define desired trajectories for the 
remaining states the tracking of which automatically ensures the original tracking objective for xi(*). Indeed, the 
intuition behind such a definition of desired trajectories becomes clear when we look at X 2 d(t). 

X 2d{t) = — - *1 [Xi (t) - X ld (t)] 

From (166) it is clear that when the surface x 2 = X 2 d the resulting dynamics for xi(f) is given as 
xi(t) = x 2 (t) 

— X2 d(t) 

= ~ M*l(*) - *ld(0] 

The dynamics of the system is such as to ensure that Zi(t) — xi d (<) exponentially. However, if the surface x 2 -z 2d = 0 
can only be reached exponentially, then the dynamics of x\ is perturbed by an exponentially decaying signal, and 
therefore invoking the result on the exponentially stable systems perturbed by exponentially decaying perturbations, 
we conclude exponential convergence of Xi(f) to x\ ,d(t). We now show the relationship between control laws developed 
using the recursively defined desired trajectories and the standard pole-placement control law. 


(166) 

(167) 

(168) 

(169) 


Theorem 8.3 Connection between pole-placement and recursive trajectory definition. 
Given (Gl) A nonlinear system of the form (161) - (164). 

(G2) Given a set of desired trajectories of the form (165) 

If (II) The controls u in equation (I 64 ) are chosen to be 


« = - M*n(0 - Xnd{t)] 


(170) 


where x,, d (t) : 5?" x 9? + -*• R »' = 2, 3, . . . , n is specified by (165) and k„ € 97+ . 

Then (Tl) The control law specified by (170) is a stable pole-placement control with the n eigenvalues each being equal 
to —k t i = 1, 2, . . . , n 


Proof: 4 t> The proof is obvious by writing the dynamics for x\ and i 2 - Indeed, 

- k 2 [x 2 {t) - x 2 d(t)\ 


X\ — X2 

dx 2d (t) 

X 2 ~ dt 

Using the definition of x 2 d(t) provided by (165), we rewrite (172) as 


(171) 

(172) 


x i = x 2 

, _ -*!-(«)]] 

Z 2 = 

-k 2 [x2(t) - [^1 (<) - Ild(l)]]] 

Which may be rewritten as 


(173) 

(174) 

(175) 


Xi = X2 

d 2 

X2 


+ [fci + + [*ifca]*id(«) 

-[*1 + ^2]x 2 (<) - [fcl* 2 ]zj(t) 


That is to say 


Xi 


0 

X2 


. + [*1 + ^ 2 ]^^ + [* 1 Jfc 2 ]x ]d (t) 


0 1 

~[£i k 2 ] — [&i + k 2 ] 



x 1 


X2 


(176) 

(177) 

(178) 

(179) 

(180) 


The placement of poles through recursive trajectory definition is trivially obvious by inspection of equation (180). 
This concludes the proof of the theorem. <!<► 
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Comment 8.6 It is to be noted that this tracking control law is valid for any specification of desired trajectories that are 
smooth, the tracking of which guarantees achievement of the control objective. That is, we are free to specify any smooth 

set of trajectories Xid(t) i = 2,3 n, the only constraint being x,(t) = x,d(t) i = 2,3, ...,n => x,-i (<) -+ *i-i, d(t). 

Indeed, the linear pole-placement control law is just a special case of control laws that achieve this tracking objective. 

Comment 8.7 We now ask if it is possible to relax the smoothness assumption on the desired trajectories x, d {t). 
Indeed, the first relaxation would be to consider desired trajectories that are differentiable almost everywhere, except 
possibly on sets of zero measure. The Nested and Switching control laws presented in the previous chapter are examples 
of such discontinuous control laws, the discontinuities existing on sets of zero measure. The proofs of such control 
laws are much harder in general, though the regularization of such control laws that involve saturation functions have 
been used in the recent literature. We have been inspired by the attempts of [Tee92] in developing control laws that use 
Filippov averaging instead of regularization. That is to say, that we are prepared to tolerate chatter and limit cycling by 
using discontinuous control laws. The drawback however is that we can show finite time synchronous stabilization only 
on the average, whereas a regularized control law, by eliminating the discontinuity would permit smooth stabilization, 
though exponentially , without the chatter. 

Comment 8.8 Our interest »n relaxing the smoothness assumption on the desired trajectories merely enables us to 
utilize the discontinuous , synchronous control law for a practical mechanical system. 


We will first present the control law for a group of n £ Z+ mechanical systems, and then apply it to a well known 
example of the magnetic bearing. Many mechanical systems are represented by Newtons force and torque balance 
equations that assume the form 


x\ = x' 2 

X 2 = V * 


(181) 

(182) 


where x % € 3R 2 is the state of the ith mechanical system where i < n £ Z+ } and v l (x,t) : 9ft 2 x 5R+ -+ is the input 
force. Typically, xi represents the generalized position coordinate of the mechanical system, and x 2 represents the 
generalized velocity coordinate. These equations, though simple in form, serve to illustrate the application ^of the 
theory, and also represent a large class of useful physical systems. Given desired trajectories x\ d (t) : 9?+ $h to be 
tracked by the states x\ (t), we attempt to find control laws u' that ensure synchronous tracking for the states £*,(<). 

We now state the theorem that ensures synchronous tracking for the systems of the form (181)-(182). 


Theorem 8.4 Synchronous tracking for a class of mechanical systems. 

Given (Gl) n mechanical systems, each of the form (181) - (182). 

(G2) Given a set of desired trajectories of the form x* ld (t) : 3R+ 5R i = 1, 2, . . . , n 
If (II) The controls u*(x, i) t = 1, 2, . . . , n in equation (182) are chosen to be 


dx o , 




x 2 


x 7d = 


dx\ d 


dt 


■K 


E:.iW-*iaP] 


Then ( Tl ) 
(T2) 


whereki , k 2 € . 

Filippov solutions exist for system (181) - (182) subject to control (183). 
States rl(t) trac ^ ^eir respective trajectories x\ d (t) synchronously. 


(183) 

(184) 


Proof: ♦> The proof is simple once we realize the validity of the system equations (181) - (182) subject to the 
control law (183) for arbitrarily small neighborhoods of the origin. Indeed, the control law is undefined only on a set 
of zero measure. As this set of zero measure is indeed the set we desire to make invariant, and the control law directs 
system trajectories to this set, and hence maintains invariance, the conclusions of the theorem naturally follow. The 
theorem can also be proved invoking the results of the nested and switching control laws mentioned in the previous 
section. 


8.3 Application to the Magnetic Bearing 

In this subsection we apply the proposed tracking control law to magnetic bearings. 
The dynamics of magnetic bearings are given by the following equations. 

i'l = A 
i' 2 = ti' 

i = 1, 2, 4 
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(185) 

(186) 
(187) 



Given desired trajectories xj d (f) i = 1, 2, ... ,4 to be tracked by the respective state variables x\ (t) 1 
We now define the following set of vectors. 

1! 

to 

e)(t) = 

x j{t) x jd{t) i — 1,2, ... ,4 j — 1,2 

(188) 

ei(i) = 

[ «!(t) e}(<) ] T 

(189) 

e 2 {t) = 

[ «2 (<) e\(t) ] T 

(190) 

where 



*2 d(t) = 

d(t) » — x\d(t) Jf || /.\j| / n 

(191) 

= 

di ^ lf ll e l(0l|2 = 0 

(192) 

Now note that x^t), * s not strictly differentiable at the origin, but has a derivative that exists almost everywhere. 
Indeed define the generalized derivative as 


dt ||e,(«)||* 

(193) 

= 


(194) 

where 



m) = 

T2 e i(<)[d(*)«i(o - «l(<)ei(<)] 1 = 1 . 2 , 3,4 j # . 

(195) 

e\(t) = 

x\(t)-x\ d {t)i = l,2,...,4 

(196) 

We now choose v 

1 i = 1 , 2, . . . , 4 in the following manner. 


v'(t) = 

±i-«) - if ||e 2 (0lb * 0 

(197) 


d*U*) 

dt 


if ||e 2 (<)l |2 = 0 


(198) 


Claim 8.1 Synchronous Tracking for a Magnetic Bearing . 

Given (Gl) Mechanical systems, each of the form (187), 

( G2) Given a set of desired trajectories of the form x\ d(t) : 9R+ — ► 3? i = 1 , 2, . . . , n 
If (II) The controls v*(x,t) i = 1, 2, ... ,n in equation (187) are specified by (197) , wherek\ , h? £ 9?+ . 
Then (Tl) Filippov solutions exist for system (187) - subject to control (197). 

(T2) States x\(t) track their respective trajectories xj d (t) synchronously. 


Proof: ^ t> The proof of the claim is by invoking the theorem proved earlier for the more general case of a group 
of mechanical systems. 

Indeed, it is easily seen that the application of control (197) would cause the states x£(<) i = 1,2, .. . ,4 to reach 
their desired values in finite time, and the desired trajectories are so chosen that the reduced dynamics ensures finite 
time tracking for xl d (f). <] + 

Results of simulation are shown for the following conditions. The chosen desired trajectories were as follows. 
r idW = sin t, x.\ d (t) = 5, Xj d (<) = —2, x} d (f) = 5. The initial conditions were as follows x{(0) = 1, xf(0) = 

7, £](0) = — 1, x}(0) = 2, 12 ( 0 ) = 0, xl(0) = 0, i 2 (0) = 0, ^(0)=0. 


Simulation results are in excellent agreement with the predicted behaviour. Indeed, note that the trajectory 
errors vanish identically at the same instant of time. This indeed was the motivation for considering the synchronous 
tracking control law. 
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Figure 6: Synchronous Control of a Magnetic Bearing 

9 Fractional Control of the Magnetic Bearing 

9.1 Introduction 

In this section, we will present an interesting variable structure control law for a vector dynamical system, that is a 
bounded control law, but whose convergence rate is faster than a comparable linear control law, and whose robustness 
properties are much better than comparable linear control laws [Pra92]. We will clarify what we mean by comparable 
linear control laws in the following subsections. We use the term fractional control law to indicate that this is a 
particular form of variable structure control law where the powers of indices are positive fractions. 

We will present qualitative arguments for the conjecture, and will provide simulation results that are in agreement 
with the conjecture. However the proof of this conjecture has been quite elusive, and we have been unable to present 
anything more tangible than this conjecture. We leave the proof of this control method as an open problem to the 
reader. 


9.2 Finite Time With Continuous Control - Scalar Systems 

Consider a scalar dynamical system of the form 

(199) 

x = u v ' 

where x € 9? and the control u € ». Given the control objective of regulating the state of the system (199) to the 
origin commencing from arbitrary initial conditions in finite time, we choose u in the following manner. 

„ = (2 ° 0) 

where k 6 $?+ and r > 1. 

Comment 9.1 The choice of u t s novel since the control is obviously continuous, but not differentiable at the origin. 
Also note that the control law involves raising the power of |x| to a fraction, and hence the term fractional control. 

We now make the following claim regarding existence of trajectories, stability and convergence for the system 
(199). 

Claim 9.1 Existence of solutions, stability and convergence for fractional control of scalar systems. 

Given 

(Gl) System dynamics of the form (199) 

If 

(II) The control u is specified as in (200) 

Then 

(Tl) Cauchy solutions exist for (199) subject to (200). 
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(T2) x = 0 *j stable. 


(T3) Indeed x — ► 0 in finite time t *, given by t* = ij 


v = ?- 

2 


Proof: 4* t> Existence of Cauchy solutions is easily seen by the fact that therighthand side of the differential 
system is continuous. 

Considering the candidate Lyapunov function V(x) : 5? — ► given by 

„2 

( 201 ) 

Indeed V = — < 0. Attractivity of the origin is therefore confirmed. 

To show finite time convergence we solve the equation 

X = -k\x\r sgn[x] (202) 

to obtain that t* = . The proof of the claim is complete. <14|l 

We now make a comparison between three kinds of control laws that regulate the state of the system (199) to the 
origin. 


^linear — 

— kx 

(203) 

sltdtng — 

-o* if I 1 ! > 0 
1*1 

(204) 

W fractional ~ 

k 

r* X if Ixl > 0 

kl 1 "* 

(205) 


Comparison of control efforts reveals something interesting. For all |x| > 1, the linear controller has the maximum 
gain, closely followed by the fractional controller, and the sliding mode controller has the smallest gain. However the 
situation is reversed when |x| < 1. . 

Similarly, the times taken to reach the origin from initial conditions x(0) ^ 0 are 


tlinear — 

gliding = 

OO 

l*(«)| 

k 

(206) 

(207) 

t fractional = 

kWI 1 -* 

(208) 

*[!-£] 


We now formulate an alternative control law that combines the best of both the linear and the fractional control 
law to give 

u* = —kx if | x j > 1 (209) 


-t— r-x if 0 < |x| < 1 
\x\P 


p > 1 


(210) 

( 211 ) 


Note that we do not bother to define the control law at the origin. 

There is yet another viewpoint as to why this control law does better than a linear control law when |x| < 1. The 
linear control law has an eigenvalue — k , and but u* has an eigenvalue (we use the term eigenvalue very loosely 

1 * 1 ' . 

here, since strictly speaking even the term eigenvalues does not make sense in a nonlinear context) that is increasing 
to oo as |x| — ♦ 0. Though both control laws are bounded, qualitatively, the fractional control law converges much 
faster to the origin as seen in the following scalar example. 

Example 9.1 Fractional Control - Scalar Case 

Consider the simple scalar example given by the equations 


x = u 
Choose 

^linear 

U* 

k 

P 


— kx 

— kx if \x\ > 1 

~j^ x '/ 0 < M < 1 

2 

2 


( 212 ) 

(213) 

(214) 

(215) 

(216) 
(217) 


It is clear from the simulation plots that the modified fractional control law outperforms the linear control law. 
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Now consider a linear system in the controllable canonical form, given by the following equations. 


X\ = X 2 


(218) 



(219) 

x n = u 


(220) 

where x € u £ 3?. 

Now choose the control u to be of the following form 


u — — k\X\ k 2 x 2 **• k n x n 

if ||z|| 2 > 1 

(221) 

k\ k 2 

- j,*l ~ *2 ’ 

• ■ k \zi if 0< ||s|| 3 < 1 

(222) 

IMl! Il*ll 2 r 

11*11! 


where 

n 

11*112 = 


(223) 

»=i 

r > n 


(224) 

s n + k n s n 1 + ki 

is a stable Hurwitz polynomial 

(225) 


We now formulate the following conjecture. 

Conjecture 9.1 Existence of solutions, stability and convergence for fractional control of controllable linear systems. 
Given 

(Gl) System dynamics of the form (218) - (220) 

If 

(11) The control u is specified as in (221) - (225) 

Then 

(Tl) Filippov solutions exist for systems (218) - (220) subject to control (221) - (225) 

(T2) x — 0 is globally stable 

(T3) Indeed x — ♦ 0 faster than a comparable linear control law of the form uunear - —k - k 2 x 2 - • ■ ■ - k n x n 

Qualitative Proof: 

First we note that within the unit ball ( 1 1 jc | ] 2 < l)i th e control effort is bounded by 

M<X> ( 226 > 

« = 1 

So the control does not blow up at any instant of time. We have used the notion that in the nonlinear setting, 

within the unit ball, we have each eigenvalue A, i = 1, 2, . . . , n of this system being replaced by — ^ where r > n. 

IMI 2 r 

Consequently, from the way the A, t = 1,2, . .. , n combine to form the k x of the control law, the form of the control 
law is intuitively obvious. 

We find by simulation that the robustness, and rate of convergence of the proposed nonlinear law are much 
superior to a linear control law. The proof of this conjecture, however, has eluded us. 


Example 9.2 Fractional Control for Magnetic Bearing 
We present simulation results for a system of the form 

*i = *2 


(227) 

(228) 


where 

V — —k\x\ - *2^2 *7 IWU > 1 

= — * 4 -*i - if 0 < iMia < 1 

11*11! 11*11# 

it, = 6 

Jt 2 = 11 


(229) 

(230) 

(231) 

(232) 
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The results show the faster convergence of the state subject to fractional control. 



Position - Linear Control 



Position - Fractional Control 



Figure 7: Fractional Control of a Single Axis of a Magnetic Bearing 


10 Sliding Mode Observers for the Magnetic Bearing 

10.1 Introduction 

We first present the basic theory of sliding mode observers for mechanical systems, and prove the existence of 
generalized Filippov solutions and stability. We then show the convergence of the observer state errors to zero. We 
then present the problem with existing theory, and present bounds on variables that would prevent observer failure. 
Finally we remark on the utilization of the computed bounds as a design rule to help design such sliding mode 
observers. 

The problem of designing observers using sliding mode theory was first introduced and studied by [Mis88]. Here 
the observation problem is treated as a special case of a state regulation problem. Sliding surfaces are designed based 
on the error dynamics, and reaching a sliding surface is equivalent to the error in the estimate of the measured state 
decaying to zero. In sliding mode control, the surface S = 0 is reached in finite time, and on that surface the states 
decay exponentially. Similarly, in sliding mode observer theory, the error in the estimate of the measured state decays 
in finite time. All other state errors decay exponentially. 

Consider a simple mechanical system of the form 


x\ = X 2 (233) 

i 2 = u (234) 

where x 6 and u £ 5?. Now consider an observer of the following form. 

x\ = X 2 4- kisgn[xi] (235) 

x 2 — k2sgn[x2] (236) 

x = x — x (237) 

Such an observer structure equation leads to error dynamics of the form 

x i = %2 — k\ sgn(x\ ) (238) 

X 2 = — ^25^n[xi] (239) 


Theorem 10.1 Convergence of the state estimation errors: 
Given 

(Gt) Error dynamics of the form (238)- (239) 

If 

(II) |x 2 | < *1 

Then 
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(Tl) Generalized Filippov Solutions exist for the system (238)- (239) 

(T2) The one-dimensional manifold x\ = 0 is attractive 

(T3) The averaged dynamics of x 2 about the surface X\ = 0 decays exponentially . 

Proof: 4 t> Existence of Filippov solutions is due to the fact that the governing differential inclusions are closed, 
bounded, convex and uppersemicontinuous. 

We will prove the theorem using simple Lyapunov analysis. Consider the candidate Lyapunov function, 



Differentiating V along the flow of the system(238), we get, 

V - xi[x 2 + kisgn[xi]] 

< -||£i||[*i " X2sgn[ii]] 


(240) 


(241) 

(242) 


Thus as long x 2 < fci, V < 0* indeed the surface x = 0 is attractive. 

Comment 10.1 The Theorem asserts the existence of a tubular neighbourhood around the x 2 = 0 axis where the 
trajectories converge to the manifold given by x i = 0. It is to be noted that x 2 must not be greater than k\ until the 
trajectories converge to x i =0. Some additional conditions are necessary to prevent such an occur? ence. 


The dynamics of the system when constrained to evolve on the surface x\ - 0, can be derived using 

the Fillipov 

solution concept. Thus, taking a convex combination of the dynamics on either side of the sliding surface, we get, 

x\ — 

i[i 2 + * 1 ] + (1 - “r)[*i - fc i] 

(243) 

X 2 = 

7*2 + (1 - 7)( — * 2 ] 

(244) 

From the above equations, we eliminate 7, and from the invariance of the sliding surface, we get, 


i\ = 

0 

(245) 

x 2 = 

*2 - 
Jfcl 12 

(246) 

Exponential decay of x 2 is clear from the above equation. The proof of the theorem is complete. <4 
We will now utilize this design technique to design sliding mode observers for the magnetic bearing, 
magnetic bearing system represented by the following equations. 

Consider the 


- 1 

*1 = 

1 

*2 

(247) 

• 1 

x 2 = 

1 

u 

(248) 

■ 2 

*1 = 

2 

X 2 

(249) 

x\ = 

u 2 

(250) 

i\ = 

A 

(251) 

■ 3 

X 2 = 

4 , 3 

— a * u> * x 2 4- « 

(252) 

• 4 

*1 = 

4 

*2 

(253) 

• 4 

X 2 = 

4 1 1 

a * u * x 2 + u 

(254) 

We design 

a sliding mode observer for this system of the following form, 


■ 1 

Xj = 

x 2 + ki sgn[x\] 

(255) 

* 1 

x 2 = 

u 1 T k 2 sgn[x\] 

(256) 

; 2 

x x = 

A 4- kisgn[xl] 

(257) 

; 2 

x 2 = 

u 2 4 - k 2 sgn[x]] 

(258) 

;3 

Xj - 

A 4- kisgn[A ] 

(259) 

;3 

X 2 = 

— a * w * x 2 4 - u 3 4 - fc 2 «^n[xi] 

(260) 

; 4 

*1 = 

x\ + fcis$rn[x}] 

(261) 

l 4 

X 2 = 

a * u) * A 4" u* 4" 

(262) 

where x* = x* 

-x) i = 1,2,..., 4 j = 1,2, and k u k 2 > 0. 
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We write the observer error equations as follows. 


1 1 
*1 

= 

*2 - fcisf/nfi}] 

(263) 

i, 1 
*2 

= 

-k 2 sgn[x\\ 

(264) 

L 2 
*1 

— 

il - *i sgn[x\] 

(265) 

L 2 

x 2 

= 

- k 2 sgn[x \ ] 

(266) 

X. 3 

*1 

= 

^2 - k 1 sgn[xl] 

(267) 

*2 

= 

—a * w * x 2 — k 2 Sgn\x\\ 

(268) 

L 4 

*i 

= 

x\ — k\sgn[x\] 

(269) 

x\ 

= 

a * u> * x\ — k 2 sgn[x\] 

(270) 


We now state the result concerning the stability and convergence of the observer states to their true values. 


Theorem 10.2 Convergence of the state estimation errors: 

Given 

(Gl) Error dynamics of the form (263)- (270) 

If 

(it) i*ii < 

Then 

(Tl) Generalized Filippov Solutions exist for the system (238)- (239). 

(T2) The one- dimensional manifold x\ =0 is attractive . 

(T3) The averaged dynamics of x 2 about the surface x\ = 0 decays exponentially. 

Proof: ^ > Stability of the error dynamics is easily shown utilizing the following Lyapunov function. 


X>l*i 

1=1 

(271) 

4 

-kik 2 y^[sgn[r)]] 2 
1 = 1 

(272) 

0 

(273) 


Furthermore, V — 0 — = 0. Invoking the invariance principle of LaSalle, it is seen that the largest invariant 

set containing the set x[ = 0 * = 1.2. ... ,4 is the set £\ = 0, * = 1, 2, . . . , 4. Stability, and hence convergence to the 
origin is therefore assured. 

We show convergence of the states x\,t = 1,2,..., 4 to the origin in finite time as follows. As the system is 
asymptotically stable, there exists an instant of time t* such that for all t > t ,taT , ||i|| 2 < Jfc,, _ |f'| < jfcj. Invoking 
the^ theorem on the finite time convergence of the sliding mode observer, the observer states converge in finite time. 


11 Closure 

We have shown a variety of nonlinear controllers for the magnetic bearing that are simple and robust to build and 
are guaranteed to be stable. We believe that the design of controllers utilizing principles of nonlinear analysis provides 
new richness, insight and excitement in the design of high precision magnetic bearings. 


References 

[AB90] Charalambos D. Aliprantis and Owen Burkinshaw. Principles 0} Real Analysis. Academic Press Inc 

1990. 

[AC84] Jean- Pierre Aubin and Arrigo Cellina. Differential Inclusions - Set- Valued Maps and Viability Theory. 
Springer- Verlag, 1984. 

[AEB75] Y. C.Ho A. E. Bryson. Applied Optimal Control. Hemisphere Publishing Corporation, 1975. 


51 



[AKP911 R. Gurumoorthy A. K. Pradeep. On concluding stability for Lyapunov analysis of max-functions. Pro- 
ceedings of VARSCON 91, Variable Structure Control of Power Conversion Systems, Reno Nevada, June 


[AKP92a] 

[AKP92b] 

[AKP 5] 

[Bar76] 

[Fil61] 

[Fil88] 


1991. 

S. M. Shahruz A. K. Pradeep. A control law that achieves simultaneous tracking for a group of systems 
in finite time. International Journal of Control, 56:1419-1428, 1992. 

S. M. Shahruz A. K. Pradeep. A control law that achieves simultaneous tracking for a group of systems 
in finite time. Proceedings of the 1992 American Control Conference, pages 3134-3135, June 1992. 

S. M. Shahruz A. K. Pradeep. Design of control laws that achieve simultaneous tracking for a group of 
systems in finite time. ESRC Technical Report - VC Berkeley, pages 3134-3135, ESRC 91-5. 

Robert G. Bartle. The Elements Of Real Analysis. John Wiley & Sons., 1976. 

A. F. Filippov. Application of theory of differential equations with discontinuous righthand sides to 
nonlinear problems. Proceedings Of The 1st IFAC Congress, 1:923-927, 1961. 

A. F. Filippov. Differential Equations with Discontinuous Righthand Sides. Kluwer Academic Publishers, 
1988. 


[KF70] 

[Kor68] 

[Mis88] 

[Pra92] 

[RG93] 


[Rud64] 

[SMS93] 

[SMS 7] 


A. N. Kolmogorov and S. V. Fomin. Introductory Real Analysis. Dover Publications Inc, 1970. 

Korivin. Planar dynamical systems. J. Math. Mech, 8:47-60, 1968. 

A. Misawa. Observation of Nonlinear Dynamical Systems. PhD thesis, Massachussetts Institute of 
Technology, 1988. 

A. K. Pradeep. Sliding Mode Control of Perturbed Nonlinear Systems. PhD thesis, University of Cali- 
fornia, Berkeley, 1992. 

A. K. Pradeep R. Gurumoorthy. Modelling and control of a rotor supported by magnetic bearings. 
Proceedings of the Second International Symposium on Magnetic Suspension Technology, Seattle WA, 
August 1993. 

Walter Rudin. Principles Of Mathematical Analysis. McGraw-Hill Book Company, 1964. 

A. K. Pradeep S. M. Shahruz. A high precision synchronizing control system for biaxial positioning 
tables. Proceedings of the 1993 American Control Conference, pages 2239-2240, June 1993. 

A. K. Pradeep S. M. Shahruz. A high precision synchronizing control system for biaxial positioning 
tables. ESRC Technical Report - UC Berkeley, March 1992, ESRC 92-7. 


rSMSarl A. K. Pradeep S. M. Shahruz. A high precision synchronizing control system for biaxial positioning 
tables. ASME Journal of Dynamic Systems, Measurement, and Control, December 1993 (to appear). 

[SRS91I B Sheng S R. Sanders, R. Gurumoorthy. Analysis technique for practical implementation of variable 
structure control systems. Proceedings of VARSCON 91, Variable Structure Control of Power Conversion 
Systems, Reno Nevada, June 1991. 

[SS831 J. J. Slotine and S. S. Sastry. Tracking control of nonlinear systems using sliding surfaces with applica- 

tions to robot manipulators. International Journal Of Control, 38. No. 2:465-492, 1983. 

[Tee92] A. Teel. Global stabilization and restricted tracking for multiple integrators with bounded controls. 

Systems and Control Letters, 18, 1992. 

[YCBDB82] C. D. Morette Y. C. Bruhat and M. Dillard-Bleick. Analysis, Manifolds And Physics. Elsevier Science 
Publishers, 1982. 


12 Appendix - Mathematical Foundations of Discontinuous Con- 
trol 


In this section we present basic definitions, facts, and examples about measures of sets and functions integrabdity, 
absolute continuity, convexity, and differential inclusions. For further details about the definitions to follow, refer to 
[KF70], [AB90], [YCBDB82], [Rud64], [Bar76], [AC84], We will use the following concepts to develop solutions ol 

differential equations with discontinuous right hand sides. 


12.1 Measure of Sets 

We commence by formalizing the notion of an interval. 

Definition 12.1 Let Si p denote p dimensional Euclidean space. By an interval in we refer to the set of points 
x — [x*j , . . . , x p ] T such that 

, . (274) 

3a, < x, < (« = 
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The possibility that a, _ 6, is not ruled out, and the empty set is also included as a possible candidate for the 
intervals. An interval can be understood to refer to a p dimensional cube in 5? p . 


Definition 12.2 If A is the union of a finite number of intervals, then A is said to be an elementary set. 


Definition 12.3 If I is an interval in 9£ p , the measure p of the interval I is defined to be 
p:ftP^SR = X p =1 (6,-a0 


(275) 


The measure of a set is in some rough sense the volume of the geometric object formed by that set. Indeed, the 
measure of a k < n — 1 dimensional object in n dimensional space is 0. Therefore, the measure of a point in is 0, 
and the measure of a plane in is also 0. The kinds of control we work with will vanish on an n — 1 dimensional 
subspace of n dimensional space. Therefore they vanish on a subspace of 0 measure. 

Definition 12.4 If the set A is the union of a finite number of pairwise disjoint intervals, (i.e) A — (J n _j h where /_, p| I k = 
0 V j ^ k, then the measure p of the set A is J 

n 

= ^ **(•*>) (276) 

j=i 

Fact 12.4>1 Open sets are measurable. 

Fact 12 . 4.2 The union of a sequence of measurable sets is also measurable. 

Fact 12.4.3 The complement of a measurable set is also measurable. 

Fact 12.4 4 A set consisting of one point is measurable. Its measure is 0 . 

Fact 12.4.5 A denumerable set (a union of a sequence of countably many one-point sets) is measurable. Its measure 
is 0 . 

Fact 12 . 4.6 Every subset of a set of measure 0 is measurable. Its measure is 0 . 


Example 12.1 Measure of Set of Undefined Control 
Consider the control u(t) : 9R+ — ► given by 

u(t) = — sgn[x] (277) 

where x € 3?. Note that the function spn[z] : SR - {0} — [-1, 1] is not defined at 0. Using fact (12.4-4), we assert that 
the control is not defined on a set of zero measure. 


Definition 12.5 By almost everywhere, we mean everywhere excepting possibly on a set of measure 0. 

Simply, by saying a relation holds true almost everywhere, we assert that the set of points where the relation fails to 
hold, has measure zero. 

Example 12.2 Behaviour of Functions 

Consider functions f : X — », g : X — f n : X — » n = 1 , 2, . . .. The following are the instances of almost 
everywhere relations between the functions. 

I f = 9 almost everywhere , if p{x £ X : f(x) ^ </(x)} = 0. 

2- / > 9 almost everywhere , if p{x £ X : f(x) < ^(x)} = 0. 

5. fn * f almost everywhere, if p{x £ X : f n (x) ^ f{x)} ~ 0. 

4- /nt/ almost everywhere, if f n < / n+1 almost erert/where for all n and 

5. fn l f almost everywhere, if / n +i < f n almost everywhere for all n and 

Example 12.3 Controls Defined Almost Everywhere 
Consider the control u(i) : — ► 5? given by 

u(i) = -s<?n[z] ( 2 78) 

where x £ 5? Note that the control u(t) is defined almost everywhere. 


fn * f almost everywhere, 
fn —*■ f almost everywhere. 
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12.2 Measurable Functions 

In the definitions to follow, we will consider real valued functions that map a measurable space A' with measure p to 
the extended real line 9ft. 

Definition 12.6 A function / : X — * 5? is measurable if the set S = {x € X : |/(j:)| < a} is measurable for all 
a > 0 € 9?- 

Measurability of a function is a property of the function, based on the measurability of a certain set in its domain. 
Example 12.4 Measurable Functions 

All continuous functions f : 9ft n — ► 9ft are measurable. Proof is by showing that by definition , S is open , and hence 
measurable. 

The following are some facts based on operations between measurable functions. 

Fact 12.6.7 Given f : X — * 9ft, if f is measurable, then Vk G 9ft, kf is measurable. 

Fact 12.6.8 Given f : X — ♦ 3?, and g : X — * 3?, if f and g are measurable, then f + g, f - g, and fg are measurable. 

Fact 12.6.9 Given p : X — », and f : X — », if f is measurable, and if f(x) = p(x) almost everywhere, then p is 

also measurable. 

We now present a fact concerning the properties of the limit function, based on the properties of the convergent 
functions. 

Fact 12.6.10 Given a sequence of functions f, : X — 3? i = 1, 2, . . . convergent almost everywhere to the function 
f : X -> if each /•,« = 1,2, .. . is measurable then the function f is also measurable. 


12.3 Integrable Functions 

The key idea of the integral developed by Lesbegue is as follows. In Riemann integration, if / : [a, b] 6 9ft — * 9ft, then 
to form the Riemann integral we divide the domain of /, [a, b] into many subintervals and group together neighboring 
points in the domain of the function /. On the other hand, the Lesbegue integral is formed by grouping together 
points of the domain where the function / takes neighboring values in the range ! Indeed, the key idea is to partition 
the range of a function rather than the domain. It is immediately obvious that such a technique allows us to consider 
functions that have multiple points of discontinuity, or may not even be defined at some points. 

The functions in this section map a measurable space A" with measure p into the extended real line 9ft. 

Definition 12.7 Let f : X -► 9ft be a measurable function . Furthermore, let f take no more than countably many 
distinct values jji , y 2 , . . . , • . * in its range. Then, the Lesbegue integral of the function f over the set A, denoted by 

f A f(x)dp. is given by 

/(.)*. <.<*,) fm) 


where 

A, = {x : x e A, f(x) = yj} ( 28 °) 

if the series (280) is convergent. If the Lesbegue integral of the function f exists, then we say the function is 
integrable, or summable with respect to the measure p on the set A. 


Example 12.5 Integrable Functions 

Consider the constant function f(x) = 1. Let us evaluate the Lesbegue integral of the function. Indeed, 

( 281 ) 
( 282 ) 

As the Lesbegue integral exists, the function is said to be integrable. 

Fact 12.7.11 Given a measurable function f : X - *, if there exists a sequence f„ : X —> of integrable functions 
converging uniformly to f on the set A, then the function f is said to be integrable. 


J f(x)dp = J 1 dp 

= P(A) 
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Fact 12.7.12 Given f : X — ► 5ft, g : X — > 3?, if g > 0, and 15 integrable on a set A, and |/(jc)| < g(z) almost 
everywhere on A, then f is also integrable on A, and 

I / < / 4>(x)dp (283) 

J A J A 

Fact 12.7.13 Given f : X — ► 9?, t/ / is bounded and measurable on a set A, then f is integrable on A. 

Fact 12.7.14 Given f : X — ► 5ft, if f is integrable on a set A, then f is integrable on every measurable subset of A. 

Fact 12.7.15 Given f n : X — ► 5ft » = 1, 2, . . a sequence of functions converging to a limit f : X — ► 5ft almost every- 
where on a set A, if there exists a function g : A' — 5ft integrable on the set A, and\f(x)\ < g(z), V x almost everywhere in A, 
then f is integrable on A, and 


lim 

n— *00 


/ fn(x)dft = / 

J A J A 




(284) 


Fact 12.7.16 Given f„ : X — ► 5R, a sequence of functions converging to a limit f : X — ► 5R almost everywhere on a set 
A, if 3A: € SR such that |/(z)| < k, V z almost everywhere in A n = 1,2,..., then f is integrable on A, and 


lim 


f fn(x)dn = f 
J A j A 


f(x)dn 


(285) 


12.4 Absolute Continuity 

Definition 12.8 The function f : X -* 5? is said to be absolutely continuous on the interval [o, 6] 6 5? if for any 
e > 0, there exists b > 0 such that for finitely many disjoint open intervals (a, , fc, ) C [«, 6] 

n n 

^ l fr - - Oi| < 6 => ^|/(fr0 - /(a,)| < e (286) 

1=1 I 

That is, the function is of bounded variation. It is to be noted, however, that functions of bounded variation need not 
necessarily be absolutely continuous. 

Fact 12.8.17 An absolutely continuous function f : X — ► 5? is necessarily continuous. 

Note however, that the converse is not true. We will now present an example to illustrate that the converse is not 
true. 


Example 12.6 Continuity and Absolute Continuity 

Consider the function f : [0, l] — + 5ft defined as follows 


/( 0 ) = 0 

f(x) = X 2 COS -ir 0<£<1 
x z 


(287) 

(288) 


The function f is differentiable at each x £ [0, 1], but is not of bounded variation. This is trivially shown by considering 
a partitioning as follows. 


Pn 


Variation of f 



oo 


(289) 

(290) 

(291) 


The /unction is not of bounded variation, and hence is not absolutely continuous. We have thus shown an example 
of a continuous function that is not absolutely continuous. 


Fact 12.8.18 If f : X — ► 5ft, then if f is differentiable, the f is absolutely continuous. 
Fact 12.8.19 Any function that satisfies the Lipschitz condition is absolutely continuous. 


The importance of absolute continuity is that for Lesbegue integration, the fundamental theorem of calculus holds 
precisely only for absolutely continuous functions. 
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12.5 Convexity 

Definition 12.9 A set M € 5R" is said to be convex if it contains the joining of any two points in the set. 

Example 12.7 Set of Bounded Functions 

Consider C[ a ^ j, the space of all continuous real valued functions f : [a, 6] 6 9£ — 1 - fef M be the subset of 

C(a,6] defined as follows 

M = {/ 6 Cj a , 6 , : |/(*)| < 1} (292) 

Then , the set M is convex . This is easily evidenced by considering the join of two functions /(x),g(x) € M For any 
a, /? > 0 such that a + fi = 1, we have 

\af(x) + pg(z)\ < |« + P\ ( 293 ) 

< 1 (294) 

6 M ( 295 ) 

The join of two elements of M belongs to M , and convexity is shown . 

Fact 12.9.20 If a set M is convex, so is its interior. 

Fact 12.9.21 The intersection of a finite number of convex sets is also convex. That is, 

Mi is convex (296) 

i=i 

where each M, is convex. 

Definition 12.10 The minimal convex set containing a convex set M is called the convex hull of M . 

12.6 Set Valued Maps 

Definition 12.11 A function f(x) : 5? n — ft is said to be upper semi-continuous at a point x* € if 

f(x*) = lim sup f(x) (297) 

Example 12.8 Upper- Semicontinuous sgn Function 

Consider the function sgn[x) : 5R — * [— 1, l] defined as follows 

sgn[x] = 1 x > 0 
= -1 x < 0 

It is obvious that at x = 0, fhe function is upper-semicontinuous. 

Definition 12.12 Having defined some properties of real-valued functions , we now move on to discussing some 
important properties of set- valued functions. Given two sets X andY , we define a map F :X to be set- valued, 

if F : X Y associates to any x 6 X , a subset F(x) € Y . 

Definition 12.13 The domain of a set valued map F : X —> Y is defined to be Domain(F) = {x eX:F(*)/0) 
and the range Range(F) of a set valued map F : X — ► Y is defined to be (J j6X F(x). 

Definition 12.14 The domain Domain(F) of a set valued map F : X — ► Y is defined to be strict if Domain(F) — 
X, and is defined to be proper if Domain(F) ^ 0. 

Definition 12.15 A set-valued map F : X -> Y is defined to be compact if its range Range(F) is a compact subset 

of Y . 

Definition 12.16 A set valued map F : X -> Y is upper semicontinuous at x* € X , if for any open set Sy 
containing F(x*) there exists a neighbourhood Sx of x* such that F(Sx) C Sy . 

In this section we present basic results for the local existence of solutions of differential equations with discontin- 
uous righthand sides. We define a sliding mode, and present conditions for the existence of a sliding mode. We then 
present briefly the development of the sliding mode control law, and the various regularizations of it. 

We will now state without proof the following two important results from analysis that we will need. 

Arzela-Ascoli Theorem: .... 

Let K be a compact subset of R p and let F be a collection of functions which are continuous on A and have values 

in R 1 . The following properties are equivalent. 


( 298 ) 

( 299 ) 
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1. The family F is uniformly bounded and equicontinuous on K . 

2. Every sequence from F has a subsequence which is uniformly convergent on K . 

The theorem allows us to define a sequence of approximate solutions of a differential equation, and guarantees 
convergence of the approximate solutions to a limit function of the sequence is equicontinuous and uniformly bounded. 

Filippov Convergence Lemma: 

Given a differential inclusion of the form x = T{x,i). If the inclusion F(x,t) is closed, bounded, convex, and 
uppersemicontinuous, the limit of any uniformly convergent sequence of approximate solutions of the differential 
inclusion ts also a solution of this inclusion, in the domain of convergence. 

That the limit function satisfies the differential inclusion is the main reason for invoking the lemma. 
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